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Abstract 

We apply methods of dynamical systems to study the behaviour of universe dominated by the 
generalized Chaplygin gas. We reduce the dynamics to a 2-dimensional Hamiltonian system and 
study its behaviour for various ranges of parameters. The dynamics is studied on the phase plane 
by using methods of qualitative analysis of differential equations. The behaviour of trajectories at 
infinity is studied in some convenient coordinates introduced on the phase plane. Hence we shown 
that FRW model with the generalized Chaplygin gas is structurally stable. We clearly find the 
domains of cosmic acceleration as well as conditions for which the horizon problem is solved. We 
also define some general class of fluids which generalize the Chaplygin gas. The dynamics of such 
models in terms of energy conditions is also discussed. 
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I. INTRODUCTION 

The generalized Chaplygin gas jlj was recently proposed as a alternative to the cosmolog- 
ical constant in explaining the accelerating universe. The equation of state, describing the 
pressure p of this gas, is given in terms of the energy density p by the relation p = —A/p a , 
where A is a positive constant and < a ^ 1. An universe filled by matter in the form of 
Chaplygin gas evolves from the phase dominated by dust to a de Sitter phase where p ~ —p, 
p = A through an intermediate regime described by the equation of state for the Zeldovich 
stiff matter with p = p. 

There are many various ways in which one can explain the transition from a matter 
dominated universe to one with the accelerated expansion, for example brane inspired models 
[3] or vacuum metamorphosis ^|. The generalized Chaplygin gas is a generalization of 
Chaplygin gas j^| motivated by considering a d-brane in a spacetime of d+2 dimensions. 
We should note that there is no such motivation for the generalized Chaplygin gas. At a 
toy model, one can simply generalize equation of state for the Chaplygin gas p = —A/p to 
the more general dependence on density p = —A/p a n]]. Interesting future of this model 
is that it is naturally provided explanation of transition from a decelerating phase at early 
stage of evolution to accelerated expansion at later stages. 

There are monting evidences that our universe accelerates at present Q. Therefore it 
should be dominated by some kind of matter with negative pressure, the so-called dark 
energy. The different dark energy models are motivated by the recent measurements of 
high-redshift supernovae la observations, which cannot be explained by the canonical E-deS 
model. The simplest model of dark energy is the cosmological constant but in this case 
we cannot explain why we don't observe the large energy density pa = — 10 76 GeV, 
expected from particle physics which is about 10 123 times larger than the value predicted by 
the FRW model. 

Therefore, while the most obvious candidate for such energy seems to be vacuum energy, 
we are looking for alternatives which satisfies the equation of state for the quintessential 
matter p = w(z)p for simplicity, where w(z) = p/p is a quintessential parameter as a 
function of redshift z. In this context we consider the universe filled by the generalized 
Chaplygin gas as a possible candidate for dark energy. 

In investigation of full dynamics we use the dynamical systems methods to analyse the 
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global properties of the system. We visualize all evolutional paths for all possible initial 
conditions as curves (and singular solutions) in the phase plane (H,p), where H = J^(lna) 
is the Hubble function. The phase space diagrams are very useful both to analyze the 
assymptotic states, their stability and to study other interesting properties of time evolution 
of models under studies. Especially, it is interesting to study the solution of horizon problem 
and initial conditions for acceleration of our universe. 

As complementary to phase space analysis, the Hamiltonian approach is developed. In 
this approach the dynamics is reduced to the 1-dimensional Hamiltonian system. Hence, 
all possible evolutional paths are classified in the configurational space {a: a ^ 0}. The 
construction of the Hamiltonian is very advantageous and can be applied not only for quan- 
tum cosmology. It is possible to make the qualitative classification of evolutional paths by 
analyzing the characteristic curve which represents the boundary of the domain admissible 
for motion of the system. 

The next advantage of representing dynamics in terms of Hamiltonian is possibility to 
discuss the stability of critical points which is based only on the geometry of the potential 
function. 

The presented formalism gives us a natural base to discuss the property of structural 
stability of the model. We give a simple proof of structural stability of the FRW model with 
the generalized Chaplygin gas. 

The structural stability is sometimes considered as a precondition of the "physical exis- 
tence" . The existence of many drastically different mathematical models well agreeing with 
observational data (taking into account final measurement errors) seems to be fatal for the 

n 

empirical methods of modern physics [13J. Therefore, any structurally unstable models are 
probably not physically meaningful (such as some of discussed in the present work). This 
suggests - if we agree that physically realistic may be only what is structurally stable - that 
the Chaplygin gas (or other component type of possitive cosmological constant), in the real 
universe, should be different from zero. In this paper we also demonstrate the effectiveness of 
analysis of dynamics in terms of one-dimensional Hamiltonian flow to answer the question: 
How are trajectories with interesting properties distributed in the phase plane? Along which 
trajectories is the acceleration condition, a = —dV/da > satisfied? Are the trajectories, 
along which cosmological problems are solved, distributed in some typical or exceptional 
way? It is the problem of degree of generality for interesting properties. 
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Organization of the text is following. In section 2 the FRW model with the generalized 
Chaplygin gas is investigated in terms of dynamical systems methods. In section 3 we find 
that the models with Chaplygin gas are structurally stable. In section 4 we consider the 
dynamics of FRW model with the generalized Chaplygin gas as a 1-dimensional Hamiltonian 
system and prove the structural stability of this model. Section 4 gives a summary of 
results and comments on them. In section 5 we discuss how the model with the generalized 
Chaplygin gas fits the SNIa data. 

II. THE FRW MODEL WITH THE CHAPLYGIN GAS AS A TWO- 
DIMENSIONAL DYNAMICAL SYSTEM 

The main aim of qualitative analysis of differential equation is investigation of global 
dynamics in the phase plane for all possible initial conditions instead finding its explicit 
solutions. The solution of dynamical system x = f(x), f G C°° is given by a function of 
two variables, an initial condition x as well as the time t — > x(xo,t). The space of all 
states x of the system at given time is called phase space. We reduce dynamics to the 
2-dimensional phase space in which singular solutions x = are represented by critical 
points and nonsingular ones by phase curves. In this representation the phase diagrams in a 
two-dimensional phase space allow to analyze the acceleration and the horizon problem in a 
clear way. We reduce the dynamics to the two-dimensional phase space with an autonomous 
sytem of equations x = P(x, y), y = Q(x,y), where x, y are coordinates, P, Q e C°°, and a 
dot denotes the differentiation with respect to cosmological time t. 

The classification of non-degenerate critical points can be given in terms of eigenvalues 
of a linearization matrix at a critical point. The eigenvalues Ai, A2 are invariants of critical 
points, i.e., they do not change as we change the coordinates x, y. We call a critical point 
(xo,yo) non-degenerated (hyperbolic), if ReAi 7^ and ReA2 7^ 0. The particular 
trajectories of the system approach a critical point for t — > 00 or escape away from it for 
t — > —00 along a direction vector k = (k x , k y ). These directions are simply eigenvectors at 
the critical point (xo, yo). If Ai 7^ A2, then a slope of a tangent vector to trajectories as they 
are reaching this critical point is k = k x /k y given in terms of eigenvalues 

Ai - P x {x ,y ) A 2 - P x (x ,y ) 

P y (x ,yo) ' 2 P y (x ,yo) {) 
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and a prime denotes the differentiaion with respect to x or y. 

The behaviour of the system in the neighbourhood of the critical point (xo,yo) is quali- 
tatively equivalent to the behaviour of its linear part 

x = K( x o, Vo)(x - x ) + Py(x , y )(y - y ) (2) 
y = Q' x (x ,y )(x - x ) + Q' y (x , y )(y - y ). (3) 

After integration, system ©-© gives 

x - x = Re(Cie Al * + C 2 e Mt ) (4) 
y-y = Re(C 1 k 1 e Xlt + C 2 k 2 e X2t ). (5) 

These formulas (JT|)-(|5)) allow to study the behaviour of the system near its non-degenerate 
critical points and then to draw global phase portraits of a dynamical system. For com- 
pleteness of analysis the behaviour at infinity is considered. For this aim it is useful to 
compactify M 2 by adding a circle at infinity. One can do that simply by introducing the 
projective coordinates. 

The two maps cover the circle at infinity x = oo, y = oo: 

z = — , u = — ; z = 0, — oo < u < +oo (6) 

X X 

and 

1 x 

v = -, w = — ; v = 0, — oo < w < +oo. (7) 

v y 

Let us consider now the FRW model with the vanishing A-term in the following repre- 
sentation 



<7 J-f 1 

H = — = -H 2 --(p + 3p) = P(H, p), (8a) 
p=^ = -3H(p + p) = Q{H,p), (8b) 

where all symbols have their usual meanings. The first integral of (|8a|) is p = 3H 2 + 3/c/a 2 , 
where k is curvature and a is a scale factor. 

The equation of state for the generalized Chaplygin gas 

P = ~, A>0, 0<a^l (9) 
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has been assumed. In equation (|8a|). (H,p) are chosen as phase variables. Because the 
dependence p on H signifies, the bulk viscosity in general p = p(H, p) can be postulated. The 
evolution of the system is represented by trajectories in the (H, p)-space uniquely determined 
by the initial conditions (H ,p Q ). 

The first equation from the system (J8a|) can be recognized as the Raychaudhuri equation 
for the congruence of the world lines of matter in spacetime with the R-W symmetry. The 
second one is a condition of the conservation energy-momentum tensor. 

Let us introduce the following sets (SEC) & (WEC) 

W = {(H,p): p + p^O}, (10a) 
S = {(H,p): p + 3p^0}. (10b) 

System (jHaj) has the critical points of two types: 

1. static (H = 0) critical points of the system (jHJ) are determined by intersection of the 
p-axis and the boundary dS of the S set, i.e. {static critical points} = {p axis} D dS; 

2. nonstatic critical points of (jHJ) are situated on the intersection of the flat FRW 
trajectory p = 3H 2 and the boundary of W, i.e. {nonstatic critical points} = 
{trajectory k = 0} H dW . 

From the physical point of view critical points on the phase space corresponds to an 
asymptotic states of the universe. To determine its character (i.e. a type of stability) we 
consider linearization of (jHJ) at the critical point (Ho,po)- The linearization matrix of the 
system is 

8P 8P 
a_ dHi dp 

dQ dQ 
_dH' dp 

The characteristic equation det(^4 — Al) = takes the form 

A 2 - (Tr^)A + detA = 0, (12) 

where 

(P + 3p) (13) 

o 



-lipip+M 

-3(p + p), -3H±(j> + p) 



(11) 



TrA = 0, detA 



1 , . d 
'2 ip + P) °Tp 
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for static critical point and 

d 



TiA = -H, 



2 + 3 J 
dp 



detA = QH%4- 
dp 



(P+P) (14) 



o 



(p + p) 



for non-static critical points; here denotes the value of the expresion at the critical point 
(H ,p ). 

Let us assume now that 1) the phase space does not contain trajectories along which the 
energy condition p + p ^ is violated in a finite value of t. In terms of dynamical system 
this assumption is equivalent to the condition that dW is a trajectory in the phase plane 
since p = on dW . We also postulate that 2) p + p ^ only for p > p mhl . On the strength 
of the above assumptions the evolution of the universe is described by one of the trajectories 
lying between p m i n and infinity and no trajectory intersects dW in a finite value of time. 
Every model violates the strong energy condition. 

In this case there is only one static critical point of the saddle type (then eigenvalues 
are real and opposite signs) and two non-static critical points, stable node if H > and 
unstable node if H < 0. They represents deS types of evolution. 

The phase diagram of such a model is shown in Fig. [IJi. For comparison the phase 
portraits for dust filled models with the cosmological constant are also presented (Fig. 
d). In these cases 

H = -H 2 -\p+^ p = -3H P . 

Let us note that all flat and open solutions experience acceleration for the late time after 
crossing line {p = p crit } which corresponds to the energy density for the Einstein static 
universe. All models which are situated close to the flat model also accelerate. 

There are some characteristic domains in Fig. but both the model with generalized 
Chaplygin gas (Fig. ^i) and the model with positive cosmological constant A (Fig. ^i) have 
equivalent portraits in their physical domains: 

I - open models are starting from a singularity and going towards the stable deS node 

II - closed models are starting from a singularity and going to the global attractor - the 
stable deS node 

III - trajectories moving in this region correspond to the closed models starting from the 
unstable deS node towards the stable deS node after bounce. 

IV - the trajectories situated in region IV are confined by the separatrix going from the 
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FIG. 1: The phase portraits for a) the FRW model with Chaplygin gas (A = I, a = 1) and the 
dust filled FRW models with cosmological constant b) A < 0, c) A = and d) A > 0. Let us note 
the equivalence of phase portraits from figures a) and d) in the domains p ^ and p 

respectively. 



singularity to the saddle point and by the separatrix coming out from the saddle point. 
This region is covered by trajectories describing the closed models evolving from the initial 
singularity to the maximal size and recolapsing to the singularity. 

V - the trajectories of region V decscribe the open models which begin their evolution 
from the unstable deS node and end on line p = on which the transformation of time is 
singular. 

VI - in this region trajectories go out of the stage p = to the stable node. 
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Let us note that because the phase portrait is symmetric when H changes its sign, in the 
regions V-VI, I- VIII and II- VII the picture is quite symmetric. 
We consider our system (jSJ) in the physical region defined as 

{{H, p): p> A > 0, < a < 1}. 

For the system under consideration the relation p(a) can be obtained from (|8b|) . Integra- 
tion gives 

Relation (|15|) interpolates between a universe dominated by dust and de Sitter one via a 
phase described by the soft matter equation of state p = jp. It would be useful to rewrite 
()15|) in the form 

p{a) = Po,cha P i I A s + fl 3 (1+a ) 1 , (16a) 

or 

pip) = Po,cha P i(A s + (1 - A s )(l + z) 3 ^)) 1 ^, (16b) 

where l + z = a~\ A s = A/(Po,ch ap i) 1+a and p ,cha P i = (A + Bfl^ +a \ 

Let us notice that even though equation admits a wider range of positive a the sound 
velocity cj. = dp/ dp = aA/p 1+a does not exceed the velocity of light. Therefore we assume 
0<a^l,£>>0 and the case of a = 1 corresponds to the Chaplygin gas. 
The first integral of (JSJ) takes the form 

3k ( 1 — A \ 3k 

3Ij2 = P ~ ^2 = P0,Chapl [ A s + Q 3 (1+a ' J - — , (17) 



a 2 



or 



h 2 = H 2 {n Chap i, (A s + (i - a s )(i + zf 1+a ^ 1+a + n kj0 (i + zf], (18) 

where f^chapi.o = Po,cha P i/3i?o ^ s the P resen t density parameter for the generalized Chaplygin 
gas. If we substitute into ([THJl z = (H = H ) then we obtain the constraint ^cha P i,o + ^fc,o = 
1. 
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Let us note that trajectories of the FRW model with generalized Chaplygin gas lie on the 
algebraic curves given by the relation 



3 M^-^J =P"3^ 2 , (19) 

where p > 3H 2 for k = +1, p = 3H 2 for k = and p < 3H 2 for k = —1; where p ^ A T +^ . 

Therefore the phase space is divided by the parabolic curve p = 3H 2 on different domains 
with respect to the curvature index. The flat model trajectory {k = 0} separates the regions 
of the model with the negative and positive curvature. 

After the substitution of the form of equation of state into (jHJ) and the reparametriza- 
tion of time t such that a new time variable is a monotonic function of original cosmological 
time we obtain the equivalent dynamical system 

d JL = _#2 pQ _ 1 ( a+1 _ 3A) = p {Hj j (2Qa) 

dr 6 

^ = -3H(p a+1 -A) = Q(H,p), (20b) 
dr 

where t — > r: dtj p a = dr. Of course relation ()19|) plays role of the first integral for (j2TJj) . 

System (|20j) is a special case of the systems for which the weak energy condition is satisfied 
only for p > p min = A l K l+a \ Therefore we have critical points of two types 

1. static H = 0, p = (3A)^ 

2. nonstatic H = ^A 5 ^), p = A^ 

The critical points lie in the physical region {(H, p): p ^ A T +^ , H e E}. From the physical 
point of view the static critical point is representing the static Einstein universe whereas 



non-static one is deS solution. For the static critical point TiA = and detA = — 



A(l+a) 

Therefore the characteristic equation for eigenvalues of the linearization matrix admits for 
this case only real solutions of opposite signs. This means that static critical point is a 
saddle point for any value of a and a positive value of constant A. For the non-static critical 
point (H ^ 0) we have: 



TrA = -H (5 + 3a), 
detA = 6# 2 (1 + a), 
A = (TrA) 2 - UetA = H 2 (9a 2 + 6a + 1) > 0. 



10 



Therefore for any a we have stable (Hq > 0) and unstable (Ho < 0) nodes. The structure 
of phase space does not depend on the particular choice of a 6 (0, 1]. 




FIG. 2: Analysis of behaviour of trajectories in phase plane transformed on the compact Poincare 
sphere for a) the FRW model with the Chaplygin gas (A = l,a = 1) and the dust filled FRW 
models with the cosmological constant b) A < 0, c) A = and d) A > 0. 

The detailed analysis of the system under consideration can be performed on the phase 
portrait (Fig. [1^). The phase trajectory for the flat model (k = 0) is given by a parabole, 
displayed symmetrically with respect to the p-axis, whose vertex changes its position de- 
pending on the value of A. The behaviour of trajectories at infinity can be investigated by 
using the mapping of the phase plane onto the "Poincare sphere" which transforms the orig- 
inal phase plane into the projective (compact) plane. If we consider the FRW model with 
the positive cosmological constant and dust then there are three critical points at infinity, 
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FIG. 3: Analysis of behaviour of trajectories for the system (|2U|) in phase plane (H,E). The 
transformation p — > E = ^fp splits the composite critical point (H = 0, p = +00) into two simple 
non-degenerate points. 



in the physical region of the phase plane (i.e. p > 0). These are two saddles (H = ±00, 
p = 0) and the degenerate point (H = 0, p = +00). The transformation E = ^fp splits the 
composite critical point into two simple non-degenerate points. 

In the case of dynamical system (j201) with the right-hand sides in the form of polynomials 
of degree n (maximal), one uses the projective coordinates. In this coordinates, points which 
at infinity correspond to a circle S l are covered by two lines z = 0, —00 < u < +00, and 
w = 0, —00 < v < +00. The dynamical system in the projective coordinates (z, u) and after 
the reparametrization of time da/dr = H n ~ l , assumes the form 

u = — = z n Q[-,-) -uz n P(-,-), 22 
da \z z J \z z J 

where n = 2 in the considered case. 

The analysis of this system is carried out in standard way jjj. The results of such an 
analysis, for dynamical system (|2"Uj) are presented on Fig. and Fig. El 
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III. STRUCTURAL STABILITY OF THE FRW MODEL WITH THE GENERAL- 
IZED CHAPLYGIN GAS 

The idea of structural stability originated with Andronov and Pontryagin p. Let us 
consider the autonomous dynamical system S' (i.e. its r.h.s) which is a small perturbation 
of original one S. The dynamical system S' is said to be structurally stable if dynamical 
systems in the space of all dynamical systems it is close, in the metric sense, to S or it is 
topologically equivalent to S. 

In the qualitative theory of dynamical systems instead of finding and analyzing an indi- 
vidual solution of a model, a space of all possible solutions is investigated. A property is 
believed to be realistic if it can be attributed to large subsets of models within a space of 
all possible solutions or if it possesses a certain stability, i.e., if it is shared by a slightly 
perturbed model. There is a wide opinion among specialists that realistic models should 
be structurally stable, or even stronger, that everything should possess a kind of structural 
stability. What does the structural stability mean in physics? The problem is in principle 
open in more than 2-dimensional case where according to Smale there are large subsets 
generic of structurally unstable systems in the space of all dynamical systems [3J. For 2- 
dimensional dynamical systems as in the considered case the Peixoto theorem states that 
structurally stable dynamical systems on compact manifolds form open and dense subsets 
in the space of all dynamical systems on the plane. Therefore, it is reasonable to require 
the model of a real 2-dimensional problem to be structurally stable. When we consider the 
dynamics of 2-dimensional models then there is a simple test of structural stability, namely 
if the right-hand sides of the dynamical systems are in polynomial form the global phase 
portraits is structurally stable S* 2 (R with a Poincare sphere) if and only if 1) a number of 



critical points and limit eye 
connecting saddle points 



es is finite, 2) each is hyperbolic and there are no trajectories 
14j . In the considered case the points at infinity are revealed on 
the projective plane. Two projective maps (z,u), (v,w) cover a circle at infinity given by 
{z = 0, — oo < u < +00} and {w = 0, —00 < v < +00}. 

The structural stability is sometimes considered as a precondition of the "real existence" . 
To have many drastically different mathematical models all of them equally well agreeing 
with the observational data (taking into account final measurement errors) seems to be fatal 

n 

for the empirical method of modern science |6(. Therefore, any 2-dimensional structurally 
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unstable models are probably not physically meaningful. These tendencies in contemporary 
cosmology inspired G. F. R. Ellis to formulate what is called by him the probability principle, 
which states that the universe model should be one that is a probable model within the set 
of all universe models and the stability assumption, closely related to the above principle, 
stating that the universe should be stable to perturbations. 

From the analysis of phase portraits at a finite domain as well as at infinity one can con- 
clude that model under consideration is structurally stable in its physical domain. Because 
there are no separatrices connecting saddle points, all critical points are hyperbolic and all 
critical points are non-degenerated and their number is finite. It should be treated as an 
adequate model of the real universe. Moreover, it is true for the more general class of models 
for which p + p ^ only for p > p min > 0. 

Let us consider now the general class of such models. Then on the strength of our 
assumption that the phase space does not contain trajectories along which the Lorentz 
invariability condition p + p(p) ^ is broken in a finite time, the phase portrait can be 
simply obtained. 

The universe evolution is described by one of the trajectories lying in the domain p > p m j n 
and there is no trajectories which intersect dS in a finite value of time. Then Universe 
accelerates because the strong energy condition is violated. Other energy conditions may be 
satisfied at all times. 

Let us note that structure of the phase space does not depend on the present form of 
pressure because the critical points and their stability depend on the energy condition. Of 
course, the obtained phase portrait is equivalent to that obtained for the Chaplygin gas 
model. 

IV. FRW MODELS WITH THE CHAPLYGIN GAS AS A HAMILTON! AN SYS- 
TEMS 

It is well known that the first integral FRW equation p — 3H 2 = 3k /a 2 can be used to 
construct a Hamiltonian function for the system. We take advantage of this feature in the 
considered model. The right-hand side of the Raychaudhuri equation ()8aj) can be expressed 
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H 



FIG. 4: The phase portrait for the general class of FRW models with condition that p + p ^ 
for p > p m i n > 0. The generalized Chaplygin gas can be treated as a special case of such fluids. 
The special form of equation of state is not required because the energy condition determines the 
critical points as well as its character. 



in terms of the scale factor a(t) as 

B J A 



a 



| a -2 2 |a 3 ( 1+a ) - 1 (Aa 3 ( 1+a > + B) ~ 1+a . (23) 



Therefore the universe is accelerating provided that a > a CI = (^) 3(1+a) . 

Equation (|23*j) can be rewritten in the form analogous to the Newton equation of motion 
with the 1-dimensional configuration space {a: a e M + } 

dV(a) 



da 

where the potential function 



(24) 



V(a) = -\pa 2 + Vb, (25) 

u 

where p(a) is given by (|15|): in any case p corresponds to the effective energy density and 
Vq = const. 

The first integral of ((211) is 

V(a) + y = V - \. (26) 
15 



Now, we construct the Hamiltonian function 

H=^ + V(a). (27) 

Therefore the trajectories of the system lie on the energy level 7i = E = const. 

It is useful to consider dynamics on the zero-energy level Ji — E = 0. Then we have 
Vq = k/2 and the Hamiltonian constraint coincides with the form of the first integral. 

Finally we obtain the potential of the system in the following form 

V(o) =--(a+ ^a 2 + - = -^a 2 . (28) 

v > 6 ^ a 3(i+a) y 2 6 y ' 

There are many advantages of representing dynamics as a 1-dimensional Hamiltonian system, 
namely 

1. we can discuss the existence and stability of critical points only in terms of geometry 
of the potential function; 

2. the presented formalism gives us a natural base to discuss (in terms of the potential 
function) how cosmological problems like the horizon problem or flatness problem are 
solved by models; 

3. due to the existence of Hamiltonian constraint the classification of all possible evolu- 
tional paths can be performed by the consideration of the equation for the boundary 
of domain admissible for motion 

V = {a: V{a) < 0}. 

Let us comment now some of advantages. 

In our case the Hamiltonian function takes the simplest form for natural mechanical 
systems (i.e. the kinetic energy is quadratic in velocities and the potential function depends 
only on generalized coordinates). Then, the possible critical points in a finite domain of phase 
space are only centres or saddles. In this case the characteristic equation is A 2 + d 2 V/da 2 \o = 
(at the critical point a = ao we have dV/da\ a=ao = and a| a=ao = 0). Therefore if 
d 2 V/da 2 \o > then critical points are centres, i.e. a model is structurally unstable. One can 
easily check that in our case V"(a) < which gives that only saddles points are admitted 
and the model is structurally stable. 
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The horizon problem is solved when a ^ const, (may be equal to zero) as we go to 
singularity a — ► 0. Then the distance to the horizon J a~ l dt diverges which means that there 
are no causally disconnected regions as a — > 0. Now, from the Hamiltonian constraint we 
obtain that V(0) = or in the vicinity of singularity V > —C 2 /2. Therefore the divergence 
of the potential function near the singularity indicates the presence of the particle horizon 
in the past. 

From the physical point of view it is interesting along which trajectories the acceleration 
condition a > is satisfied. One can easily observe this phenomenon from the geometry of 
the potential function. If the potential function is a decreasing function of the scale factor 
then the universe accelerates. Therefore the accelerating region is situated to the right from 
the saddle point or for the Lemaitre-Einstein universes the acceleration begins in the middle 
of quasistatic stage. In order the solution of the flatness problem means that at large a the 
universe accelerates, such that we can find constant C that a ^ C (then a > at large a). 
Hence V < —C 2 /2 at large a from the Hamiltonian constraint. 

Let us note that near the initial singularity V(a) oc 1/a like for the dust matter and V 
goes to minus infinity as we are going to the initial singularity. Hence, the horizon problem is 
not solved by the model. On the other hand for large a, V(a) oc a 2 and the flatness problem 
is solved. Let us now concentrate on the possibility of classification of evolutional paths by 
analyzing the characteristic curve which represents the boundary equation of the domain 
admissible for motion of the system in the configuration space. Because a 2 = — 2V(a) the 
trajectories of the system lie in the region D = {a: V(a) ^ 0}. The boundary of this domain 
is dV = {a: V(a) = 0}. 

From (|28J) the constant A can be expressed as a function of a as 



The plot of A (a) for different a and B is shown in Fig. El Finally we consider the 
evolutional paths as levels of A = const. > and then we classify all evolutions modulo 
their quantitative properties. Let us note that Hamiltonian (J27|) can be rewritten in the new 
form using dimensionless quantities 




B 



(29) 



x = — , t —> T = \H \t. 
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FIG. 5: The plot of A{a) (formula (|29jl) for the classification of models with potential function (|2Sj) 
for a) B = 1, k = +1, a = 0.2,0,7,0.9,1, b) a = 1, k = ±1, B = 1,1.5,3,10, c) B = 1, k = 0, 
a = 0.1,0.4,0.7,1 and c) a = 1, k = 0, B = 1,2.5,5,10. By considering levels of constant A(a) 
we obtain qualitative evolution paths. The shaded region is nonphysical because the Hamiltonian 
constraint. 

Then the basic dynamical equations are 

y = + ^cha P i,o Us + \ x 2 = -V(x), (30) 

Of course, the above dynamical system has the Hamiltonian 

n = ^ + v(x), 

which should be considered on the zero-energy level. 

As an example of application of these equations, to study the problem of cosmic accel- 
eration, consider the case of f2 fcj0 and fichapi.o ^ 0. It emerges that at present our universe 
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accelerates provided that x > for x = 1, i.e. A s > 1/3. In general the domain of accelera- 
tion depends only on the A s , namely 

-(1 - A s )x~ 3{1+a) +2A S > 0. 

Therefore there is the minimum value of x for which the universe accelerates, i.e. 

1 - A a 



3(l+a) 



X > 



2A 



This value can be always expressed in terms of redshift z 

2A S 



z < 



I- A 




V. THE REDSHIFT-MAGNITUDE RELATION FOR THE MODEL WITH GEN- 
ERALIZED CHAPLYGIN GAS 

The cosmic distance measures like the luminosity distance sensitively depend on the 
spatial geometry and the dynamics. Therefore, the luminosity depends on the present density 
parameters of different components and their form of equation of state. For this reason a 
redshift-magnitude relation may be used to determine the best parameters for the FRW 
model with generalized Chaplygin gas. 

Let us consider an observer located at r = at the moment t = to which receives a 
light ray emitted at t — ti from a source of the absolute luminosity L located at the radial 
distance n. The redshift z of the source is related to the scale factor at the two moments 
of evolution by 1 + z = a(to)/a{t\) = a /a. If the apparent luminosity of the source as 
measured by the observer is I then the luminosity distance di of the source is defined by the 
relation 

<=4 < 32 > 

where di = (1 + z)aoTi. The luminosity distance in the flat FRW model with generalized 
Chaplygin gas takes the form 

d (z) = 1 + Z [ Z = (33) 
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From above expression we obtain the relation between the apparent magnitude m and ab- 
solute magnitude M in the form 



m — M. = 5 Ioe 



where 



'1 + z) 



dz' 



+ Z'Y + ^Chapl, (^ + (1 " A a )(l + 3)3(1+*))1/(1+*) 



A* = M-51ogi7 + 25. 
In order to compare with the supernove data, we compute the distance modulus 



jjL = 5log(d L ) 



25, 



(34) 
(35) 

(36) 



where di, is in Mps. We define effective rest-frame B- magnitude m£ ff for 54 supernovae which 
relate to the HzST results through m£ ff = M& + fio, where M& is the peak B-band absolute 
magnitude of a standard supernova. 



a) 



b) 



a 0.50- 




o.oo^ 

0.50 0.60 0.70 0.80 0.90 



1.60 0.70 0.SC 

As 



0.90 1.00 



FIG. 6: Confidence levels on the plane (a,A s ) for a) O, m ,o = 0.05 and b) f2 mi o = 0.3, (£l m ,o + 
Ochapl.o = !)• The figures show the areas of the preferred value of a and A s . The shaded areas 
show the parameter regions with confidence level 68.3% (95.4%). 
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FIG. 7: Relative magnitude (residuals) with respect to the Einstein-de Sitter model for two cases: 
the Perlmutter flat model (highest curve), and the flat model with generalized Chaplygin gas 
(middle curve) for O mj o = 0.05, fichapi,o = 0.95, a = 0.51, A s = 0.73, \ 2 = 53.66 (best fit 
parameters). For £l m fl = 0.3 we obtain a = 0.95, A s = 0.95, % 2 = 53.64. We fit the model 
with generalized Chaplygin gas to the sample C of Perlmutter SNIa data. Sample A is the entire 
data set. Sample B excludes four outliers: SN1992bo and SN1992bp (the most significant outliers 
from the average ligt-curve width) with lower redshifts and two with higher redshifts, SN1994H 
and SN19970 (the largest residuals from \ 2 fitting). Sample C further excludes two very likely 
reddened supernovae, SN1996cg and SN1996cn. 




The best fit parameters are obtained by minimizing the relation 

^ = ^ K,-j,.l (37) 

In above expression, is the measured value, f/ Qi is the value calculated in the model 
described above, a^ i is the dispersion in the distance modulus due to peculiar velocities of 
galaxies 3] • We test the model with generalized Chaplygin gas using sample C of Perlmutter 
SNIa data. We may use the method of maximum likelihood parameter estimation on this 
data to estimate the cosmological parameters of interest, namely pair (a,A s ). The result 
of statistical analysis are presented on Fig. El Figures illustrate the confidence levels as a 
function of a and A s for a) Q m ,o — 0.05 and b) Q m .o = 0.3, for the flat model (Qk,o — 0). 

In Fig. [7| we present the plot of residuals of redshift-magnitude relationship for the su- 
pernovae data. 
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FIG. 8: Levels of constant x 2 on the plane (a,A s ) (M = —3.39) with a) Q m fi = 0.05 and b) 
^m,o = 0.3 for the flat model (tlk,o = 0). The figures show the preferred pairs (a,A s ). 

In Fig. |S] we show levels of constant x 2 on the plane (a, A s ), for a) f2 m o = 0.05 and 
b) ^ m ,o = 0.3, calculated as the lowest value of x 2 f° r each pair of values (a, A s ) for 
(M. = —3.39). These figures show the favoured regions of pairs (a,A s ). 



VI. CONCLUSIONS 



In the present paper we reduce the dynamics of an universe filled with the Chaplygin gas 
to the form of a particle of unit mass in one-dimensional potential V(a), where a is the scale 
factor. Such a procedure gives at once an insight into the possible evolutional paths, similarly 
as in the classical mechanics. We consider the complementary description of dynamics of the 
model on 2-dimensional phase plane. From the theory of qualitative analysis of differential 
equations we obtain the visualization of the system evolution in the phase plane (H, p) , 
where H is the Hubble function and p - energy density of the generalized Chaplygin fluid. 
Such a geometrization of dynamics is very useful to analyze the asymptotic states and their 
stability and to study other interesting properties of dynamical systems. Especially, it is 
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interesting to study the solution of the horizon problem and initial conditions for present 
acceleration of our universe. We have a neat interpretation of a domain of acceleration as 
a domain in configuration space where the potential function decreases. Therefore from 
the observation of the potential function, we can find the acceleration domain a > a cr with 
a cr : dV/da\ aci = 0. On the other hand the horizon problem is solved if V(a) is bounded as 
a — > 0. If we find trajectories for which V(a) — > — oo as a — >• then the particle horizon in 
the past is present in the model. We demonstrate that the FRW model with the Chaplygin 
gas does not solve the horizon problem whereas the flatness problem is here resolved. 

There is an opinion quite widely spread among specialists that physically realistic models 
of the world should be structurally stable or, even stronger, that everything that exists should 
possess some kind of structural stability . The universe certainly exists but its stability 
properties are by no means clear. However In the case of 2-dimensional systems there is a 
simple test of structural stability and moreover from the Peixoto theorem such systems form 
open and dense subsets in the space of all dynamical systems on the plane, i.e. structurally 
stable systems are generic whereas structurally unstable systems are exceptional (or non- 
generic). If we consider FRW models with cosmological constant then only models with 
positive cosmological constant are structurally stable (Fig. [TJI and Fig. Ell). We prove that 
the FRW models with the Chaplygin gas are structurally stable in their physical domains. 
This indicate that this model, from the theoretical point of view, should be seriously treated 
as a possible candidate to describe the actual universe which accelerates. We also discuss 
how the FRW model with the generalized Chaplygin gas fits the SNIa data. 
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